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Abstract 

We study the impact of nuclear effects on the extraction of the weak-mixing angle sin^ dw from 
deep inelastic (anti-)neutrino-nucleus scattering, with special emphasis on the recently announced 
NuTeV Collaboration 3a deviation of sin^ 6\y from its standard model value. We have found that 
nuclear effects, which are very important in electromagnetic deep inelastic scattering (DIS), are 
quite small in weak charged current DIS. In neutral current DIS processes, which contain the weak 
mixing angle, we predict that these effects play also an important role and may dramatically affect 
the value of sin^ 0\y extracted from the experimental data. 
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I. INTRODUCTION 



It is a common belief that the standard model (SM) is a low energy remnant of some 
more fundamental theory. In fact, the recent observation of neutrino oscillations was the first 
unambiguous indication of the presence of physics beyond the SM, which allows for non-zero 
neutrino masses and lepton flavor violation processes, forbidden in the SM. This evidence 
stimulates further searches for deviations from the SM predictions in physical observables 
related to other sectors of the SM. 

The precise determination of the weak-mixing angle sin^ 6^/ plays a crucial role in test- 
ing the SM of electroweak interactions. Its present value is consistent with all the known 
electroweak observables [Q]. 

A recent announcement by the NuTeV Collaboration |0| on a 3cr deviation of the value 
of sin^ 6w measured in deep inelastic neutrino-nucleus (Fe) scattering with respect to the 
fit of the SM predictions to other electroweak measurements may be the sign of new 
physics beyond the SM. This result takes into account various sources of systematic errors. 
However, there still remains the question of whether the reported deviation can be accounted 
for by SM effects not properly implemented in the analysis of the experimental data. In the 
present paper we examine the role of nuclear effects on the extraction of sin^6'vi/, since an 
iron nuclear target was actually used in the NuTeV experiment. 

The observables measured in this experiment are ratios of neutral (NC) to charged (CC) 
current events, related by a sophisticated Monte Carlo simulation to siv^Oy^r. In order to 
examine the possible impact of nuclear corrections on the extraction of sin^6'vi/, we study 
the corresponding ratios 

of neutral current (NC) to charged current (CC) neutrino (anti-neutrino) cross sections for 
a nuclear target A. As is known, neglecting nuclear effects for an isoscalar target, one can 
extract the weak-mixing angle by using the Llewellyn-Smith relation [^: 
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written in terms of NC and CC (anti-)neutrino-nucleon cross sections. Here, 



^ cos2^vfM|' a{iyf, + N ^ fi- + X) 2' ^ ' 

However, actual targets such as the iron target of the NuTeV experiment, are not always 
isoscalar, having a significant neutron excess. In addition, nuclear effects including the EMC 
effect, nuclear shadowing and Fermi motion corrections are known to be very important for 
electromagnetic structure functions. These nuclear effects may also modify the CC and 
NC structure functions, and therefore a detailed study of these effects on the extraction of 
the weak-mixing angle is essential. In principle any nuclear model which can successfully 
explain the EMC effect in deep inelastic muon-nucleon scattering can serve for our purposes. 
For definiteness, in the present work we use a particular nuclear re-scaling model We 
make complete estimates of nuclear effects on the ratios R^^'^'^ and on the resulting values of 
siv? 6w In order to reduce the uncertainties related to sea quarks, Paschos-Wolfenstein |^ 
suggested to extract sin^ 6^ from the relationship 



aiu, + N^^^-+X)-aiV, + N^^^+ + X)-P ^2 " ^^J ' 

Inspired by the Paschos-Wolfenstein relation, we will also examine nuclear effects on sin^ 9w 
by the following observable for the scattering off a nuclear target A, 



^ a{p^ + A^li- + X)-a{V^ + A^li+ + Xy 

Below we present a detailed analysis of nuclear effects starting with a brief summary of 
the formalism we will use. 



II. NUCLEON STRUCTURE FUNCTIONS 



In the quark-parton model the nucleon structure functions are determined in terms of 
the quark u{x, Q^), d{x, Q^), s{x, Q^), c(x, Q^), b{x, Q"^) and gluon g{x, Q^) distribution func- 
tions, which satisfy the QCD Q^-evolution equations. Below we only collect expressions 
for the relevant proton structure functions. The corresponding neutron structure func- 
tions can be obtained from the proton ones by the replacements u{x,Q^) < — > d{x,Q'^), 
u{x, Q"^) < — ^ d{x, Q"^). 
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The structure functions (SF) of CC reactions v{v)N I {l^)X are given by 

fY^^' = u{x){\V^,\' + \Vus\')+umVub\'e{x,-x) (7) 
+ d{x)\V^d\' + d{Q\Vj^e{x,-x) 
+ s{x)\Vus\^ + s{^,)\V,sMx,-x) 
+ c(x)(|K,p + \V,s\^) + c{^b)\V,t\Mx, - x) 
+ b{x)\Vub\' + b{a\V,b\'e{x,-x), 

here Vij are the Cabibbo-Kobayashi-Maskawa mixing matrix elements. The variable 

^(x(i+§y {k=c,b), 

\ X, [k = u, d, s), 

and the step functions 6{xc — x), 6{xh — x) take into account rescaling due to heavy quark 
production thresholds. 

The structure functions F^^^ and F^^^ are obtained from by the replacements of 
the quark distribution functions q{x, Q^) indicated in the curly brackets: 

(x, Q') = F^^ix, Q'){q{x, Q') ^ 2xg(x, Q^), Q') - 2^9(6, Q')}, (8) 
fY"'{x,Q') = 2 Fr"(x,g2){g(x,g2) - -g(x,Q2)}. (9) 

The structure functions of the NC reactions v{v)N v{v)X are 



= \{\{glf + {g\f]{u{x)+u{x) + c{x) + c{x)) (10) 

d\2 I t d\2 



+[{9vy + {girmx) + d{x) + s{x) + ^(x))}, 



F^" = 2xFf^ (11) 

Fs^' = ^[g^glHx) - u{x) + c{x) - c{x)) (12) 
+ gUl^d^x) - d{x) + s(a;) - s(x))]. 
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In the SM the vector and axial-vector quark couphngs are given by 

^ ~ 2 ~ 3 ^ ~ ~2 3 ~ 2' ~ ~2' 

In our analysis we adopt the CTEQ6 Set-2 parton distribution functions of Ref. [Q. 

III. NUCLEAR EFFECT ON NUCLEON STRUCTURE FUNCTIONS 

Here, we summarize our approach for calculating the structure functions of a given nuclear 
target starting from the free nucleon quark distribution functions discussed in the previous 
section. 

A. Nuclear Parton Distributions in the extended x— rescaling Model 

Since the discovery of the EMC effect in deep inelastic muon-nucleus scattering, various 
models P, ^ have been proposed for its explanation. In our present work, we choose one 
of the successful EMC models usually referred as the extended x— rescaling model 0]. We 
extend this model to describe the CC and NC structure functions in (anti-)neutrino-nucleus 
deep inelastic scattering. 

Let K^^j^s^{x, Q"^) = xpA(N){x, Q^), p = V, S, Ghe the momentum distributions of valence 
quarks(V), sea quarks(S) and gluon (G) in the nucleus A (or nucleon N), respectively. 

The X— rescaling model |Q is based on the fact that in a nucleus the Bjorken variable 
X of the nucleon structure functions is rescaled due to the binding energy of nucleons in 
the nuclear environment. In Ref. 0] it was found that the universal x— rescaling violates 
conservation of nuclear momentum. In the extended x— rescaling model this problem is fixed 
by introducing different x— rescaling parameters for the momentum distributions of valence 
quarks and sea quarks (gluons)in the nucleon structure function, i.e., 

K^^'\x,Q') = K]^^'\6y^s)X,Q'). (13) 

For the x— rescaling parameters of valence quarks, sea quarks and gluons we take the values 
6v = 1.026 and 6s = = 0.945 0]. Actually, only one of these two parameters is free while 
the other can be determined from the momentum conservation sum rule. 
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B. Nuclear shadowing of sea quark and gluon distributions 



The x-rescaling with the above given parameters explains the EMC effect in the region 
of medium values of x. As is well known, in the low x region the effect of nuclear shadowing 
must also be taken into account. This amounts to a depletion of the nuclear structure 
functions at low x, due to the destructive interference of diffractive channels induced by 



final state interactions In this picture, shadowing corrections to structure functions 

are the same for neutrino or for charged lepton scattering. At low we expect other 
contributions, such as those present in a vector-meson-dominance model, and in principle in 
this case there could be differences between the shadowing corrections to structure functions 
in neutrino and in charged lepton scattering. Nevertheless, a detailed analysis with all the 
important contributions taken into account has been performed [0 , with the conclusion that 



the total shadowing in neutrino induced reactions is comparable in magnitude to shadowing 
in charged lepton induced reactions. Therefore for our purposes it is enough to take the 
shadowing corrections used for nuclear structure functions in lepton charged scattering, and 
apply them to the nuclear structure functions of neutrino induced reactions. For definiteness 
we use the approach presented in Ref. in which this effect is incorporated into the 
structure functions by introducing a nuclear shadowing factor in the sea quark and gluon 
g distribution functions 

q,{x, Q2) Rf^(x) ■ q,{x, Q^), g{x, Q^) i?l(a;) ■ g{x, Q^). 

Here we assumed that the gluon and the sea quark distribution functions receive the same 
nuclear shadowing. For the nuclear shadowing factor we use the parameterization proposed 



in Ref. 12 



1 + alnAln(x/0.15), (x < 0.15), 

Rfhi^) = { 1 + 61n Aln(x/0.15) ln(x/0.3), (0.15 <x< 0.3) (14) 
1, {x > 0.3). 

This parameterization gives the summary of those features of nuclear shadowing which 
were important for explaining this nuclear effect. The parameters a, b in Eq. (0) can be 
determined from the experimental data on nuclear shadowing in ^^Fe and ^°Ca [0, |13|. We 
find the values a = 0.013 and b = —0.02. 
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C. Nuclear structure functions and Fermi motion of nucleons 



We introduce the average nucleon structure functions in a nucleus A in a conventional 

way as @] 

K^^\^,Q') = J[F^A{x,Q') - l{N-Z){Fnx,Q')~F[{x,Q'))], (15) 

where i = 1,2, 3, and A = Z + N with and Z being the number of neutron and proton in 
the nucleus A. The functions FiA{x, Q"^) are the nuclear structure functions. The second term 
compensates for the neutron excess in a nucleus A. The functions Ff{x, Q^) and -F"(a;, Q"^) 
are free proton and neutron structure functions, calculated according to Eqs. (0)-(0). The 
nuclear structure functions Fj^(x,Q^), with Fermi motion corrections in a nucleus, can be 



written as a convolution 



A 

Here the bound nucleon structure functions F^^^^ for the nucleon in the single-particle 
nucleon state with wave function ip\{p) are given by Eqs. (|7|)- ([T^) with the substitutions 
introduced in the previous section, in order to incorporate x-rescaling and nuclear shadowing. 

In Eq. (^) we use the variable z = {j)Q+pz)/f^N,Vo = + ^x, where ex is the binding 
energy of a nucleon in the single-particle state A. The single-particle wave function ip\{'p) of 
the nucleon in momentum space satisfies the light-cone normalization condition: 

"^'^ \Mv)?z=l. (17) 



(27r)3 

In the following calculation, e\ and ipxip) are taken from Ref. 0]. 

Previously, in Ref. [^, the above described approach was applied to the analysis of the 
electromagnetic structure functions of deep inelastic muon-nucleus scattering in order to 
explain the EMC effect 0. In Fig. |l], we show the results of this analysis in the form of a 
ratio of the nuclear structure functions F2 of ^^Fe and deuteron. The experimental points 
correspond to the data of the EMC Collaboration 0. As seen from Fig. |^ the theoretical 
curve derived in the adopted approach fit the experimental data with good precision. Fig. |I| 
also shows that nuclear effects on the structure functions are significant. 

In order to further check our nuclear model, we also apply it to weak charged current 
neutrino- iron DIS, which do not depend on the weak mixing angle. In this case there 
are no free parameters. The comparison with data is shown in Fig. 0. We see that the 
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FIG. 1: The comparison of our results for ratio of the nuclear structure function of ^^Fe to that 
of the deuteron with the experimental data. = 20.0 GeV^ is used. The experimental data are 
taken from Ref. 0, IS ] . 



results with nuclear corrections (solid lines) for high and medium values of x are in excellent 
agreement with experimental data [jl5l. At small x the agreement is a bit worse but still 
quite reasonable . We also show the results without nuclear corrections (dashed lines). The 
surprising conclusion is that these nuclear corrections are negligible in charged current DIS, 
even at large x ~ 0.65 values, where in the electromagnetic case there is a large effect. At 
small X the data is not very precise, but the trend indicates that the shadowing corrections 
are negligible, and there is even the possibility of antishadowing. This region certainly 
deserves further experimental and theoretical analysis. 

Having our approach verified in the case of the electromagnetic and weak charged pro- 
cesses (structure functions) we apply it to the analysis of nuclear effects in neutral current 
(anti-) neutrino- nucleus deep inelastic scattering, and study their impact on the extraction 
of the weak- mixing angle siia^ 9w 



IV. NUCLEAR EFFECTS ON THE EXTRACTION OF sin^ 9w 



In our numerical analysis we study the influence of nuclear effects on the extraction of 
sin^ 6w from the observables R'^^'^^ and i?^, taking into account some kinematical cut-offs 
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FIG. 2: The comparison of our results for the charged current DIS differential cross sections at 
E,y = 150 GeV with the experimental data [^]. Our results are given for y > 0.31 since the 
CTEQ6 parton distributions which we used are available only for > 1.3 GeV^. 



specific for the NuTeV experiment. 

The differential cross sections for CC and NC (anti-)neutrino-nucleus deep inelastic scat- 
tering, in terms of the structure functions defined in Eq. (f^), are given by [llH 



jxy^Ff *(^)(a:, Q^)+ (18) 
dxdy 71 I 



for the CC reaction, and 



^ = V £.,F{V-F'f-"(:t,0')+ (19) 



9 



for the NC reaction. 

In the event selection, the NuTeV Collaboration applied the cut off 

20GeV < Ecai < ISOGeV, (20) 

for a visible energy deposit to the calorimeter E^ai- The lower limit ensures full efficiency of 
the trigger, allows for an accurate vertex determination and reduces cosmic ray background. 

Therefore we calculate the observables It}''^ and i?^ imposing the same cut off on the 
energy Eh of the final hadronic state X, assuming Eh = Ecai- Since Eh ^ we can write 
the kinematical variables averaged over the (anti-)neutrino flux as 

x= < 1, y = _^<l. (21) 

2MNEcal 

For the average energies of the neutrino and antineutrino beams we take the values < E^ >= 
120 GeV and < E^, >= 112 GeV, as in the NuTeV experiment 



The cut off, Eq. (pO]), characteristic for the NuTeV experimental events, does not exclude 
the region of small values of where the QCD parton picture is not really applicable. 
Therefore it is not possible to calculate the total cross sections of the CC and NC reactions 
for the ratios R'^'^^ and i?^ in a theoretically controllable way. Given that we study these 
ratios at fixed values of Q^, 



= 2mN{E,(j,))xy, (22) 

which allows us to examine the dependence of the nuclear effects. This method may be 
helpful if the experimental events concentrate around some known average value of (Q^). In 
the NuTeV experiment (Q^) ~ 20GeV^, but this value was obtained from the Monte Carlo 
event simulation. The actual kinematics of the selected CC and NC events is poorly known, 
except for the above mentioned energy deposit cut off (|20D . Despite the fact that the average 
is 20 GeV^, a substantial fraction of events may correspond to relatively low values of 
Q^. Therefore our results are not directly applicable to the NuTeV experimental data, but 
only indicate some general features of nuclear effects on the extraction of sin^ 6w from i?^*^^^ 
and Rj^ relevant to this experiment, namely that in general they are sizable. The actual role 
of nuclear effects can only be revealed by their inclusion in the corresponding Monte Carlo 
event simulation. 

The results of our analysis are summarized in Tables | and ||, and graphically in Fig. |^, 
where we define 

8R = R\- Rj^, with ^ = ^/,z7, -. 
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TABLE I: Nuclear effects on R'^ and R'^ , and Jsin^ 9w extracted from R 





SR" 


6R^ 


6 sin^ 9w 


3.0 


-0.001805 


-0.000651 


-0.00301 


6.0 


-0.000646 


-0.000659 


-0.001050 


9.0 


-0.000018 


-0.001081 


-0.000018 


12.0 


0.000328 


-0.001811 


0.000600 


16.0 


0.000639 


-0.002787 


0.000960 


18.0 


0.000747 


-0.003195 


0.001140 


24.0 


0.000995 


-0.004460 


0.001488 


30.0 


0.001092 


-0.005585 


0.00164 



TABLE II: Nuclear effects on R , and S sin^ 9w extracted from R 





6R- 


6 sin^ 6w 


3.0 


-0.000878 


-0.000980 


6.0 


-0.000165 


-0.000100 


9.0 


0.000198 


0.000210 


12.0 


0.000465 


0.000509 


16.0 


0.000792 


0.000840 


18.0 


0.000912 


0.000950 


24.0 


0.001346 


0.001408 


30.0 


0.001680 


0.001780 



as differences between the NC/CC ratios calculated with (-R^) and without (-R]v) nuclear 
corrections. The quantity 6 sin^ 9w in Tables | and |I|, and in Fig. ^, are the net effect of 
nuclear corrections on siia^ 9w extracted from (solid curve) or R (dashed curve). The 
ratio makes no appreciable influence on sin^Ow, being weakly sensitive to its variations 
in the region of its possible values, and therefore the corresponding values of 6 sin^ 9w are 
not presented in Table 1. 

We estimate 6sirL^9]y in the following way. First, we calculate the NC/CC ratio R'^ 
without nuclear effects, taking for sin^ 6w the central NuTeV value sin^ 6w = sin^ 6\^^ = 
0.2277 Then we calculate the ratio i?^ with nuclear effects in the way described in 
section 3, fitting sin"^ 9^ in order to get the value of i?^ equal to R'^ calculated in the first 



11 



step without nuclear corrections. Thus we obtain the values of S sin^ 9w from the equations 



The same procedure is applied for the estimation of the nuclear correction 6 sin^ 6w from 
the Paschos-Wolfenstein ratio i?^. 

Therefore in order to get the actual value of sin^ 6w, which is the value sin^ 6\^^ extracted 
with nuclear effects, the value sin^ 6[^^ , obtained without these effects, must be shifted by 
dsin^Ow Note that the value sin^ 6^^^ corresponds to the case reported by the NuTeV 
Collaboration. In fact, the NuTeV analysis is based on the parton model equations (7)-(12) 
for neutrino-nucleon CC and NC DIS, applied to an iron target nucleus. In this approach 
one assumes that the parton distribution functions (PDF) are the effective PDF in iron 
which absorb all the nuclear target effects, and are the same for any DIS process with the 
same target. The effective PDF in iron were extracted from CC data ||l5l and then used to 
calculate the NC structure functions. However as it follows from our analysis the PDF per 
nucleon extracted from the CC data should be very close to those of a free nucleon, since 
nuclear effects in CC DIS are negligible. Thus the analysis based on on the so extracted 
PDF deals with practically free nucleon PDF risking to lose the nuclear effects which do 
not manifest themselves in the CC DIS but can be important in the NC neutrino-iron DIS 
process. Recall that these effects are very important in charged lepton DIS processes. 

From Tables |I| and 0, and Fig. ^, it is seen that nuclear effects on i?'^, R'^, R and sin^ 9w 
are large and strongly depend on the value of Q^. Moreover, at certain values of the 
nuclear correction SsirL^Ow changes its sign. These transition values are ~ 10 GeV^ 
and — 8 GeV^ for the extraction methods using R'^ and i?^ respectively. In the NuTeV 
experiment the average value of is about 20 GeV^, which lies in the region of positive 
values of 6siia^9w- For this reason it might be thought that nuclear effects enhance the 
NuTeV deviation of sin^ 9w from its SM value. However as we already noted, a substantial 
fraction of the NuTeV events may have quite small values of providing negative nuclear 
correction 5sin^ 9w to the weak-mixing angle sin^ 9w- 

The following note is also in order. Given the kinematical cut off in Eq. (^^, the Q^- 
region of negative 6sin^9\Y corresponds to the values of the Bjorken variable x < 0.25, 
where we expect strong nuclear shadowing effects. This is the case for the parameterization 
of the shadowing factor given in Eq. (|l^ . In our analysis shadowing is a dominant nuclear 
effect in the region < 5 GeV^ corresponding to x < 0.13. In more sophisticated models 



i?:;(sin2^i^))=i?5;(sin^C)), 
5 sin^ 9w = sin^ 9i^^ — sin^ 6'^ ■* 



(23) 
(24) 
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FIG. 3: The dependence of the nuclear modification to the weak mixing angle S sin^ 9w extracted 
from (solid curve) and from R~ (dashed curve). 



shadowing may be accompanied with antishadowing, acting in such a way that in the CC DIS 
their common effect is small due to self-compensation while in the NC DIS this effect could 
be larger than in our model, enhancing the nuclear corrections to the weak mixing angle 
6 sin^ 6w at the small x-values. Thus a more careful study of these effects in neutrino-nucleus 
scattering is required. 

Recently, in Ref. it was observed that corrections from higher- twist effects of nuclear 
shadowing to sin^ 9w extracted via the Paschos-Wolfenstein relation, may well be of the 
same size as the deviation from its global fit value reported by the NuTeV Collaboration M . 



V. SUMMARY AND CONCLUSIONS 



In this work, we have shown that taking into account nuclear effects such as nuclear 
x-rescaling, nuclear shadowing and Fermi motion in the nucleon structure functions of CC 
and NC (anti-)neutrino-nucleus scattering, may significantly affect the extracted value of 
the weak-mixing angle sin^ 9]y. 

The procedure used by the NuTeV Collaboration p| in order to take into account nuclear 
effects was to extract effective nuclear parton distribution functions from charged current 
DIS data, which then were used in order to obtain the value of sin^ 9]y from their data. 
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Although in principle this is a reasonable idea, in practice our results indicate that it might 
be difficult to implement, since nuclear effects play a minor role, if at all, in charged current 
DIS, while in NC DIS they are important. Unfortunately there is no direct way of comparing 
our theoretical predictions with the results of the NuTeV Collaboration presented in the 
form of ratios of NC to CC experimental event candidates with poorly identified kinematics. 
Uncertainties in kinematics are pertinent to the experiments measuring NC (anti-)neutrino 
scattering since the final state neutrinos are not detectable. As a consequence, among the 
NuTeV events there might be a substantial fraction of low events, despite the fact that 
in this range their probability decreases with decreasing . Thus in theoretical estimations 
of ratios like W ^ , R~ , the kinematical integration in the total NC and CC cross sections 
extends to the low region including = 0. This region poses the problem of calculating 
the structure functions in a theoretically controllable way, since the QCD parton model, 
presently the only ffimly motivated one, is not applicable below ~ few GeV^, where 
QCD has a strong non-perturbative behaviour. 
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